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PREFACE 


This  user's  guide  documents  CFRAG,  a  computer  program  for  the  seepage 
analysis  of  confined  groundwater  flow  of  finite  depth. 

Funding  for  the  development  of  the  program  and  preparation  of  the  user's 
guide  was  provided  to  the  Automatic  Data  Processing  (ADP)  Center,  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  Miss.,  by  the  Civil 
Works  Directorate  of  the  Office,  Chief  of  Engineers  (OCE) ,  U.  S.  Army,  under 
the  Geotechnical  Aspects  of  the  Computer-Aided  Structural  Engineering  (CASE) 
Project . 

Specifications  for  the  program  were  prepared  by  the  members  of  the  CASE 

Task  Group  on  Geotechnical  Aspects  of  CASE.  Members  of  the  task  group  during 

the  development  of  the  program  were  as  follows: 

Dr.  Roger  Brown,  South  Atlantic  Division 

Mr.  Larry  Cooley,  Vicksburg  District 

Mr.  Frank  Coopinger,  North  Atlantic  Division 

Mr.  Richard  Davidson,  OCE 

Mr.  Ed  Demsky ,  St.  Louis  District 

Mr.  Lavane  Dempsey,  St.  Paul  District 

Mr.  Phil  Napolitano,  New  Orleans  District  (Chairman) 

Mr.  Bill  Strohm,  Waterways  Experiment  Station 
Mr.  Charles  Trahan,  Lower  Miss.  Valley  Division 
Mr.  Tom  Wolff,  St.  Louis  District 

The  program  was  written  by  Mr.  Michael  E.  Pace  with  technical  assistance 
from  Mr.  Reed  L.  Mosher,  Computer-Aided  Design  (CAD)  Group,  ADP  Center,  WES. 

The  report  was  written  by  Michael  E.  Pace  under  the  guidance  of 
Mr.  Mosher  and  Mr.  Dennis  R.  Williams,  CAD  Group,  ADP  Center,  WES.  Additions 
and  comments  were  supplied  by  Mr.  Williams,  Mr.  Mosher,  and  Mr.  Thomas  Wolff, 
St.  Louis  District.  All  hand  computations  and  comparisons  were  performed  by 
Mr.  Pace.  The  work  was  managed  and  coordinated  by  Dr.  N.  Radhakrishnan, 
Chief,  ADP  Center,  WES,  and  CASE  Project  Manager.  Messrs.  Richard  Davidson 
and  Donald  M.  Dressier  were  the  OCE  points  of  contact. 

Commander  and  Director  of  WES  during  the  development  of  the  program 
and  publication  of  this  report  was  COL  Tilford  C.  Creel,  CE.  Technical 
Director  was  Mr.  Frederick  R.  Brown. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 


(metric)  units  as  follows: 

Multiply 

By 

To  Obtain 

cubic  feet  per  day 

0.02831 

cubic  meters  per  day 

feet 

0.3048 

meters 

kips  (1000  lb  force) 

4.448222 

kilonewtons 

kips  (force)  per  square  foot 

47.880263 

kilopascals 

pound-force 

4.448 

newtons 

pounds  (force)  per  square  foot 

47.880263 

pascals 

pounds  (mass)  per  cubic  foot 

16.018463 

kilograms  per  cubic  meter 

SEEPAGE  ANALYSIS  OF  CONFINED  FLOW  PROBLEMS 
BY  THE  METHOD  OF  FRAGMENTS  (CFRAG) 


PART  I:  INTRODUCTION 

Purpose  of  Program  CFRAG* 

1.  This  report  describes  CFRAG,  a  computer  program  designed  for  the 
analysis  of  confined  groundwater  flow.  Flow  computations  are  performed  using 

the  method  of  fragments,  first  introduced  by  N.  N.  Pavlovsky  in  1935**  and 

i* 

later  presented  by  M.  E.  Harr.  *  The  program  will  provide  the  following 
information: 

<i.  The  flow  per  unit  width  through  the  specified  soil  medium. 

_b.  The  head  loss  per  fragment. 

£.  The  exit  gradient  for  certain  fragment  types. 

cl.  The  resultant  uplift  and  lateral  forces  acting  on  the 
structure. 

e_.  A  sketch  of  the  structure  with  uplift  and  lateral  pressure 
diagrams  drawn  appropriately. 

Organization  of  Report 

2.  The  remainder  of  this  report  is  organized  as  follows: 

a.  Part  II  presents  an  overview  of  the  method  of  fragments  and  the 
analysis  procedure  employed  in  the  program. 

_b.  Part  III  describes  the  user’s  guide  for  data  input  and  for 
execution  of  the  program. 

c.  Appendix  A  presents  example  computer  runs  and  hand 
verifications . 


*  Three  sheets  entitled  "Program  Information"  have  been  hand-inserted  inside 
the  front  cover  of  this  report.  They  present  general  information  on  the 
program  and  describe  how  it  can  be  accessed.  If  procedures  used  to  access 
this  and  other  CORPS  library  programs  should  change,  recipients  of  this 
report  will  be  furnished  a  revised  version  of  the  "Program  Information." 

**  N.  N.  Pavlovsky.  1956.  Collected  Works,  Akad.  Nauk,  USSR,  Leningrad, 
t  M.  E.  Harr.  1962.  Groundwater  and  Seepage,  McGraw-Hill  Book  Company, 

New  York,  pp  159-165. 

t  'f  _ .  1977.  Mechanics  of  Particulate  Media,  McGraw-Hill  Book 

Company,  New  York,  pp  17 1—174 . 


Appendix  B  gives  additional  example  problems  and  presents  a 
comparison  of  solutions  using  the  method  of  fragments,  flow- 
nets,  and  a  finite  element  method  computer  program.  The 
method  of  creep  is  also  used  in  the  comparison  of  uplift 
pressures . 

Appendix  C  presents  a  simple  hydraulic  structure  and  illus¬ 
trates  the  effect  on  solutions  obtained  by  the  various 
methods  when  the  geometry  of  the  structure  is  altered. 


PART  II:  ANALYSIS  PROCEDURE 


Basic  Assumptions 

3.  The  method  o£  fragments  is  an  approximate  analytical  procedure  for 
computing  groundwater  flow.  —The  principal  assumptions  used  in  the  derivation 
of  this  method  are  as  follows: 

a.  The  flow  is  confined  and  of  finite  depth. 

lb.  Darcy's  law  is  valid^  therefore,  laminar  flow  exists. 

c_.  A  steady  state  flow  exists. 

<i.  The  soil  medium  is  homogeneous  and  isotropic. 

ei  Equipotent ial  lines  at  certain  locations  of  the  flow  region 
can  be  approximated  by  vertical  lines , ^Figure  1). 


Figure  1.  Approximation  of  actual  equipotential 
lines  by  vertical  lines 


Development  of  Equations  of  Flow  and  He ad  Loss 
for  _the  Method  of  Fragments 


4.  The  method  of  fragments  describes  the  flow  through  a  porous  medium 
by  the  expression 


where 


(1) 


q  =  quantity  of  flow 
k  =  coefficient  of  permeability 


b 


h  =  total  head  loss 
1  =  dimensionless  form  factor 

5.  In  the  analysis  of  seepage  by  flownets,  the  form  factor  : 
corresponds  to  the  ratio  N  yN .  (Figure  2),  and  the  total  flow  is  expressed 
by 

0  ■ k  ntst  - kh  r  (2 

e  f  e 

in  which 

N  =  the  total  number  of  equipotent ial  drops  occurring  within  the 
flow  region 


N,_  =  the  total  number  of  flow  channels  corresponding  to  N 


/iili 


Figure  2.  Simple  flownet  divided  inti'  fragments 

6.  Groundwater  flow  is  analyzed  using  the  method  of  fragments  by 
assuming  that  a  flow  regime  can  be  represented  as  an  assemblage  of  smaller 
flow  regimes,  or  fragments,  for  which  form  factors  are  known. 

7.  From  Figure  2  we  observe  that  the  flow  across  any  equipotent ial  line 
is  constant  and  is  equal  to  the  total  quantity  of  flow  for  the  system. 

If  the  flow  region  in  Figure  2  is  divided  along  the  equipotent ial  line 
extending  from  the  tip  of  the  sheetpile,  the  flow  region  is  separated  into 
two  fragments.  Since  the  fragment  boundaries  consist  of  equipotent iai  lines, 
the  flow  through  each  fragment  must  be  identical  and  equal  to  the  total 
quantity  of  flow  through  the  system. 


Extending  this  concept  to  a  more  general  flow  system  divided  into  n  frag¬ 
ments,  we  observe  that 


q  =  qL  =  q2  ...  =  q.  ...  =  qn  (1) 

8.  In  Figure  2,  N  /N^.  represents  a  form  factor  for  the  entire  flow 
region.  Similarly  for  each  fragment,  the  sum  of  the  equipotent ial  drops 
divided  by  the  number  of  flow  channels  in  the  fragment  is  the  form  factor  for 
that  fragment.  Therefore,  the  head  loss  in  any  fragment  may  be  expressed  as 


h.  =  f 

i  k  i 


where 


h. 

i 


i 

9. 


head  loss  in  ith  fragment 
form  factor  for  ith  fragment 


The  head  losses  for  the  fragments  shown  in  Figure  2  are, 


h,  = 


a 

k  'i 


and 


llo  = 


3.  •_ 

k  ‘2 


Solving  for  q  and  combining  with  equation  (3) 


khi  kh2 


‘  1  '2 

For  a  more  general  case  containing  n  fragments 


kh  kh  kh . 

12  l 


kh 


n 

n 


(A) 


(5) 


10.  For  the  total  t low  region  the  quantity  ol  flow  may  be  expressed  as 


h  =  total  head  loss 


<I>T  =  sum  of  the  individual  form  factors  for  the  flow  region 

11.  If  equations  (5)  and  (6)  are  combined,  the  head  loss  for  the  ith 
fragment  is 


h.  = 

l 


(7) 


This  relationship  may  also  be  related  to  flownets  as 


h.  = 

l 


(8) 


where 

N  =  equipotential  drops  within  fragment  i  . 
i 

12.  If  a  flow  region  can  be  divided  into  fragments,  for  which  the  form 
factors  are  known,  then  the  head  loss  for  each  fragment  can  be  determined. 
Once  the  head  loss  for  each  fragment  has  been  determined,  the  pressure 
distribution  and  quantity  of  flow  can  also  be  obtained. 


Fragment  Types 

13.  The  method  of  fragments  is  a  method  of  solving  the  general  flow 
equation  by  dividing  the  flow  region  into  fragments  for  which  form  factors 
are  known.  Unfortunately,  except  for  the  simplest  of  fragments,  form  factors 
are  extremely  difficult  to  determine.  However,  Pavlosky*  was  able  to  develop 
mathematical  solutions  for  a  number  of  fragment  types  by  assuming  that  the 


* 


Op.  cit.  ,  p  4. 


equipotential  lines  on  fragment  boundaries  could  be  approximated  by  vertical 
lines. 

14.  In  the  following  paragraphs  a  brief  derivation  of  six  form  factors 
is  presented.  The  content  of  this  section  is  based  on  the  method  of  fragments 
as  presented  by  Harr*. 

Fragment  type  1 

15.  Fragment  type  1  (Figure  3a)  is  a  region  describing  parallel  flow 
between  two  horizontal  impervious  boundaries.  From  Darcy's  law  for  a  flow 
system  of  unit  area  a  , 


q 


=  kia  =  k  I 


a 


which  when  rearranged 


q 


kh  r  =  kh 

L  L/a 


From  the  general  flow  equation  we  see  that  the  form  factor  is 


(9a) 


(9b) 


<t> 


(10) 


16.  For  an  elemental  section  (Figure  3b) 


d<f 


dx 

y 


(ii) 


This  elemental  section  will  be  used  to  derive  the  form  factors  for  fragments 
4 ,  5 ,  and  6 . 


*  See  Harr  (1962),  Appendix  B,  for  a  table  of  elliptic  integrals  of  the 
first  kind. 


i r  -  r  r  *.  * 


and 


m  =  sm 


TTS 

2T 


(13) 


(U) 


The  form  factor  can  be  expressed  as 

$ 


K 

K 


(15) 


Fragment  type  3 

19.  Fragment  type  3  (Figure  5)  is  an  impervious  boundary  of  length  b 
and  embedment  s  in  a  pervious  layer  of  thickness  T  .  Like  fragment 
type  2,  it  is  also  an  entrance  (Figure  5a)  or  exit  (Figure  5b)  fragment. 


rM 


f«-  b  ->j 


/"/'  /  7  /  /7  /  7  7  7  7 
b.  Exit 


Figure  5.  Fragment  type  3 


20.  The  equation  for  the  form  factor  is  the  same  as  for  fragment 
type  2.  Therefore,  the  form  factor  may  also  be  expressed  as 


4> 


(16) 


The  elliptic  integrals  K  and  K'  are  computed  using  the  modulus 


m  =  cos 


ns  J 
2T  \ 


,  2  nb  ,  2  ns 

tanh  —  +  tan 


(17) 


Fragment  type  4 

21.  Fragment  type  4  (Figure  6)  is  an  impervious  boundary  of  length  b 
and  embedment  s  in  a  pervious  layer  of  thickness  T  .  A  type  4  fragment  is 


12 


normally  considered  an  internal  fragment.  On  very  rare  occasions,  it  may  be 

used  as  an  entrance  or  exit  fragment.  However,  in  CFRAG  the  type  4  fragment 
is  an  internal  fragment  only. 


Figure  6.  Fragment  type  4 

22.  For  the  exact  solution  of  a  type  4  fragment,  the  form  factor  is 

4  ■  r  <18 

where  K  and  K'  are  elliptic  integrals  as  previously  defined,  but  with  a 
modulus  given  by 

m  =  Asn^-  A,  A^  (19 

In  this  equation 

A  =  complete  elliptic  integral  of  the  first  kind  of  modulus  A 

A*  *  complete  elliptic  integral  of  the  first  kind  of  complementary 
modulus  A' 

also 

A  T 
A'  ”  b 

The  term  sn  ^  A,  A^  is  defined  as  an  elliptic  function  of  ^  A,  A^  . 

23.  To  simplify  the  calculations  for  this  fragment,  an  approximate 
solution  proposed  by  Pavlovsky  is  used  in  CFRAG.  From  Figure  7a  the  quan¬ 
tity  of  seepage  above  streamline  AB  was  observed  by  Pavlovsky  to  be  of 
small  order  compared  to  the  total  flow.  Neglecting  this  flow,  Pavlovsky 


divided  the  flow  region  into  two  parts  labeled  active  and  passive  (Fig¬ 
ure  7b).  The  dividing  line  for  the  two  parts  is  line  CD  at  an  angle  0  . 
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Figure  7.  Flow  system  for  approximate  solution 
of  a  type  4  fragment 


24.  On  the  basis  of  his  studies,  Pavlovsky  chose  0  =  45°.  Depending  on 
the  ratio  of  b  to  s  ,  two  form  factors  are  possible. 

25.  For  the  case  of  b  _<  s  (Figure  8a)  the  active  zone  is  composed  of 
elements  of  a  type  1  fragment  of  width  dx  (Equation  11).  Therefore, 


$  =  ;b  =  /b  dx 
o  y  oa  +  x 


Integrating,  the  form  factor  is 


,{.  = 


(20) 


For  the  case  of  b  >  s  (Figure  8b) 


rs  dx  ^b  djc 

o  a  +  x  s  T 


and  by  integration 


v  =  In 


(21) 


Notice  that  in  the  case  of  b  >  s  the  first  term  of  the  equation  defines  a 
type  4  fragment  of  b  s>  while  the  second  term  defines  a  type  1  fragment. 


The  two  fragments  are  joined  at  line  DE. 


a. 


H 

f* - b 


Fragment  type  4 


Fragment  type  5 

26.  Fragment  type  5  (Figure  9)  contains  two  equal  impervious  boundaries 
of  embedment  s  in  a  pervious  layer  of  thickness  T  and  is  separated  by  a 
horizontal  impervious  length  L  .  Observe  that  fragment  type  5  is  symmetric 
and  is  twice  that  of  fragment  type  4.  Using  the  same  approximations  as  for 
the  type  4  fragment,  two  cases  are  possible:  L  £  2s  and  L  2  2s  . 


r^b— r~b-i 


Figure  9.  Fragment  type  5 


For 


L  ^  2s 


and  for 


L  >  2s 


<t>  =  2  1 


n  (X  +  f) 


L  -  2s 


For  the  case  of  L  ^  2s  ,  the  second  term  of  the  equation  provides  a  type  1 
fragment  that  is  located  between  two  type  4  fragments. 

Fragment  type  6 

27.  Fragment  type  6  (Figure  10)  is  the  same  as  fragment  type  5 
except  that  the  embedment  lengths  are  unequal.  Again,  two  cases  are  possible 
L  >  (s'  +  s")  and  L  <  (s'  +  s")  . 


ifi 


I  a" 


////////  r 


r - n 

«5'1\  S" 

r  ; 

a  L  b'4«-  6" 


a.  L  >  (s'  +  s")  b.  L  <  (s'  +  s") 

Figure  10.  Fragment  type  6 

If  the  approximations  used  in  fragment  type  4  are  applied,  the  form  factor 
for  L  2  (s'  +  s")  (Figure  10a)  is 


*  dx  .L-s"  dx  ,L  dx 

<t  =  /o  a'+x  s'  T  l-s"  a"  +  L  -  x 
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Integrating  each  term 


$  =  In 


(s’  ♦  s") 
T 


Similarly  for  L  <  (s'  +  s")  (Figure  10b) 


dx  ,L  dx 

a’  +  x  \  a"  +  L  -  x 


and  by  integration 


where 


$ 


b' 


L  +  (s'  -  s") 
2 


b" 


(24) 


(25) 

(26) 

(27) 


Pressure  Calculations 


28.  The  head  loss  for  each  fragment  can  be  determined  from  factors 
previously  developed.  From  Bernoulli's  equation,  the  pressure  head  at  the 
entrance  and  exit  boundaries  of  each  fragment  is  obtained;  and  by  assuming  a 
linear  distribution  of  the  head  loss  within  each  fragment,  the  pressure  head 
may  also  be  calculated  for  points  along  the  interior  of  the  fragment. 

29.  Recall  from  Bernoulli's  equation  for  laminar  flow  that  the  total 
head  is  defined  as 

h  =  z  +  —  (28) 

Y 

where 

z  ■  elevation  head 

P  =  pressure 

Y  =  unit  weight  of  water 
and  that  between  any  two  points 


(29) 


1  2 

z.  +  —  =  z„  +  —  +  Eh 
1  Y  2  y  L 

where 

Eh^  =  rhe  sum  of  head  losses  between  1  and  2 

30.  Based  on  equations  (28)  and  (29),  the  pressure  may  be  calculated  at 
any  point  along  the  fragment  boundaries. 

31.  In  CFRAG  an  arbitrary  datum  is  chosen  as  the  top  of  the  headwater. 

If  this  datum  is  used,  Bernoulli's  equation  reduces  to 

P2 

z,  +  —  +  Eh,  =0  (30a) 

2  y  L 

and 

P2  =  (-z2  "  IhL)Y  (30b) 

32.  In  Figure  11,  the  head  loss  gradient  R  is  approximated  by  a  linear 
distribution  along  the  boundary  ECC'E'  .  Thus, 


R  = 


m 


L  +  s'  +  s" 


(31) 


where 

h^  =  head  loss  for  the  fragment 


Figure  11.  Head  loss  within  a  typical  fragment 
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Exit  Gradient 


33.  The  exit  gradient  may  be  calculated  if  the  last  fragment  on  the 
downstream  side  of  the  defined  system  is  a  type  2  fragment.  The  exit 
gradient  as  computed  by  CFRAG  was  derived  assuming  an  infinite  depth  of 
permeable  soil.  The  exit  gradient  derived  in  this  manner  will  yield  a  con¬ 
servative  answer. 

34.  The  exit  gradient  is  defined  as 


I 


E 


h  7T 

2KTm 


(32) 


where 

h  =  head  loss  through  fragment 
K  =  elliptic  integral  of  the  first  kind 
T  =  depth  of  flow  region 
m  =  modulus 
and  the  modulus  is 

irS  .  /  ,2  nb  ^  _  2  irS 

m  =  cos  ^  W  tanh  +  tan  ^  (33) 

35.  Structural  configurations  which  do  not  have  a  type  2  fragment  on 
their  downstream  side  would  generally  be  avoided.  Without  a  type  2  fragment 
as  an  exit  fragment,  there  is  no  embedment  and  subsequently  a  lack  of  confine¬ 
ment  at  the  edge  of  the  upper  impervious  boundary.  Gradients  and  seepage 
velocities  theoretically  approach  infinity  at  this  point  of  discontinuity. 

High  local  gradients  and  lack  of  embedment  could  initiate  piping  or  "roofing" 
at  this  point. 


Applicability  of  the  Method  of  Fragments 

36.  The  method  of  fragments  is  well  suited  for  seepage  analysis  where 
the  physical  assumptions  of  the  method  (paragraph  3)  are  satisfied  and  the 
following  situations  apply: 

■a.  The  geometry  of  flow  regime  is  simple,  i.e.,  composed  of  hori¬ 
zontal  and  vertical  boundaries. 
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t>.  The  true  equipotential  lines  do  not  deviate  significantly  from 
the  fragment  boundaries. 

£.  A  reasonably  accurate  solution  is  desired  in  a  short  time  and/or 
at  a  low  cost. 


37.  This  last  situation  includes  the  following  class  of  problems: 

a.  Sizing  of  structural  features,  such  as  base  width,  location 
and  depth  of  a  sheetpile,  and  any  other  structural  features 
that  may  require  trial  solutions. 

J).  Obtaining  pressures  on  structures,  such  as  floodwalls  and 

weirs,  where  the  results  will  be  used  for  the  analysis  of  loads 
on  structures  in  conjunction  with  experience-based  factors  of 
safety. 

c^.  Determining  flow  quantities. 

As  was  stated  in  paragraph  3,  the  soil  medium  is  assumed  to  be  homogeneous 
and  isotropic.  However,  a  homogeneous,  anisotropic  system  may  be  transformed 


into  an  equivalent  homogeneous,  isotropic  system  by  applying  the  transforma¬ 
tion  of  X  =  x\]K  /K  in  the  x  direction  or  Y  =  y\lYL  /K  in  the  y  direc- 
*  y  x  '  f  x  y 

tion.  The  equivalent  coefficient  of  permeability  for  a  homogeneous,  anistropic 
section  is  K  =  ^  KxKy  '  Some  nonhomogeneous  systems  may  be  solved  by 
applying  a  method  described  by  M.  E.  Harr.* 

38.  A  practical  application  of  the  method  of  fragments  is  in  the  estima¬ 
tion  of  uplift  pressures  for  the  design  of  T-walls.  The  method  is  particularly 
useful  in  the  process  of  sizing  the  key  and  base  components  and  developing  a 
functional  relationship  between  the  uplift  pressure  and  the  base  width  or  key 
depth. 


39.  An  example  of  how  a  T-wall  might  be  modeled  is  shown  in  Figure  12. 
Since  there  is  essentially  no  flow  in  the  backfill,  the  flow  is  neglected  in 
this  area.  The  key  width  is  also  neglected  to  avoid  unnecessarily  constrain¬ 
ing  the  flow.  The  base  width  accounts  for  the  zero  key  width  because  it  is 
equal  to  the  total  width  of  the  T-wall. 

40.  The  method  of  fragments  may  also  be  used  to  calculate  exit  gradients 
by  dividing  the  head  loss  in  the  exit  fragment  by  the  depth  of  its  embedment. 
The  exit  gradient  may  also  be  calculated  by  dividing  the  head  loss  in  an 
infinitesimal  element  on  the  downstream  face  by  its  length.  The  latter 
method  is  used  in  the  program. 


*  Op.  cit.,  p  4. 
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■■ 

Figure  12.  Analysis  of  T-wall  by 
the  method  of  fragments 


Accuracy  of  the  Method  of  Fragments 


41.  The  method  of  fragments  will  yield  reasonable  results  for  applicable 
problems.  Several  factors  affect  the  accuracy  of  the  solutions.  To  ensure 
valid  results,  problems  should  be  modeled  so  that  the  flow  lines  are  not 
constricted  into  an  unnatural  path.  The  flow  through  each  fragment  type 
follows  the  general  shape  of  the  fragment.  If  several  fragments  are  used  to 
model  a  relatively  small  portion  of  a  structure,  the  flow  may  be  forced 

into  an  unnatural  path.  Modeling  of  small  areas  may  result  in  apparent  flow 
where  a  dead  area  of  virtually  no  flow  exists. 

42.  Another  factor  to  be  considered  is  the  accuracy  of  the  assumption  of 
vertical  equipotential  lines  at  the  fragment  boundaries.  Laplace's  equation 
can  be  -satisfied  both  globally  and  fragmentally  only  if  the  chosen  fragment 
boundaries  are,  in  fact,  equipotential  lines.  This  occurs  only  in  simple, 
symmetric  structures,  such  as  a  single  embedded  sheetpile.  The  greater  the 
deviation  from  the  actual  equipotential  lines,  the  greater  the  degree  of  error 
in  the  solution.  For  many  practical  problems,  however,  the  assumption  of 
equipotential  lines  coincident  with  fragment  boundaries  yields  reasonable 
results. 

43.  Solutions  also  become  less  accurate  as  the  number  of  fragments 
joined  together  increases.  This  results  from  the  inaccuracies  in  the  assump¬ 
tions  of  the  method. 


FAR.' 


i  l-l'iDl. 


S  o  u  r  c  e  .)  f  I  npu  t 


44.  Input  data  raav  be  supplied  f ri-ro  a  prepared  data  ‘  ile  or  from  the 
user's  terminal  during  execution.  If  the  d  u  a  are  input  from  the  terminal, 
the  user  may  enter  data  by  us  !•.,>  key  command  words  or  by  following  a  prompting 
sequence . 


'data  Format 


45.  Ail  input  d.t t a ,  whether  supplied  iron  a  data  file  or  from  the 


terminal , 


are  read  in  free  lie  Id  lormat. 


In  addition: 


a.  Data  items  must  be  separated  by  one  or  more  blank  spaces 
(commas  are  not  allowed  as  delimiters). 

b.  Integer  numbers  must  be  in  nondecimal  form. 

c.  Real  numbers  may  be  in  either  decimal,  nondecimal  form,  or 
F  1 o  rma  t . 

d.  User  responses  to  all  requests  for  program  control  may  be 
abbreviated  by  the  first  letter  of  the  v»ord.  For  example, 
in  response  to  the  question, 

DO  YOU  WANT  TO  Pi  or  VATKk  DKIiSSURFS?  YKS  OR  NO, 
the  user  nav  enter  Y  or  N  . 


Data  Knt rv  1  rum  Term!1  a  1 

46.  There  are  throe  ways  to  enter  data  from  the  terminal.  They  consist 
of  a  prompting  sequence,  a  conversational  mode,  and  a  nonconversat iona 1  mode. 
For  the  user  who  is  uni  ami  liar  with  t  he  program,  a  prompting  sequence  is 
available  which  vi  '■  1  uuorv  >  !•<>  use:  a!!  nvcessar-.  data ,  An  explanation 

of  the  prompting  sequence  is  ;  in  par. .graph  ”4 .  '-or  the  more  experienced 

user ,  data  may  be  e:  -<i  •> .  ::  r.J  w  •  loll  cw  oc  D  the  accompanying  data. 

When  using  command  w.-rtis  mer  ,u:t  either  a  conversational  or  a 

nonconversat  iona  !  mode  o:  data  ent  rv.  iue  >nve  >  sat  i .  a  -t  1  mode  is  for  the 
user  who  is  i. .miliar  wit  h  t  hr  command  v.-t  <!  .  bu'  .  ,mnot  remember  the  variable 
list  associated  with  them.  !hi  ,  met  bed  is  explained  in  paragraph  66.  The 
nonconversat  iona  !  r-.sie  is  :  r  ,.»•  expo;  vd  user  and  is  explained  in 
paragraph  65. 


47.  There  are  several  options  which  may  be  requested  to  provide  varia¬ 
tions  on  terminal  input.  The  options  are  explained  in  paragraphs  70 
through  77. 

48.  For  a  list  of  the  command  words  see  paragraph  73. 

Data  Entry  from  File 

49.  Data  may  be  entered  from  a  prepared  data  file.  The  data  file 
format  is  the  same  as  the  format  for  the  nonconversational  mode.  The  user 
simply  types  in  a  command  word  and  the  corresponding  data.  All  lines  of 
input  for  the  data  file  must  begin  with  a  line  number. 

50.  As  many  problems  as  desired  may  be  contained  in  a  single  data 
file.  The  end  of  data  input  for  each  problem  is  identified  by  the  command 
word  END.  After  one  file  is  exhausted,  the  user  is  given  the  chance  to 
enter  additional  runs  from  either  the  terminal  or  from  another  data  file. 

51.  The  command  word  EXIT  may  also  be  used  to  terminate  a  run.  However, 
if  at  any  time  the  command  word  EXIT  is  entered,  either  from  a  data  file  or 
from  the  terminal,  program  execution  will  terminate. 

Input  Description 


General 

52.  The  following  is  an  explanation  of  the  command  words,  requirements, 
and  variables  for  data  input.  If  data  are  input  from  a  data  file,  then  all 
lines  of  input  must  be  preceded  by  a  line  number. 

53.  The  input  information  is  divided  into  the  following  sections: 

a.  Title  of  run. 

b.  Units. 

c.  Water  data. 

d.  Fragment  data. 

e.  Termination. 

54.  Input  data  may  be  entered  in  any  order;  however,  for  some  of  the 
available  options,  certain  command  words  must  be  entered  before  the  fragment 
data.  These  options  are  explained  in  paragraphs  70  through  77. 
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55.  Any  unit  may  he  specified,  but  the  same  unit  must  be  used  throughout 
the  problem.  The  default  system  ol  units  is  in  feet,  seconds,  and  pounds. 

56.  In  the  following  description,  [ LN ]  is  used  to  denote  the  need  for  a 
line  number  when  data  are  entered  from  a  predetermined  data  file.  Single 
quotes  ('NN')  denote  the  use  ot  alphanumeric  information;  underscore  denotes 
the  minimal  amount  of  characters  required  (FRAGMENT) .  All  command  words  may 
be  abbreviated  to  the  first  four  characters-. 

Input  information 

5 7 .  Title  of  run  (opti ona 1 ) . 

a.  [LN]  1  NAME '  -  tit  in  (h()  characters  or  less) 

b.  Note:  If  sever.:!  problems  are  run  and  a  new  title  is  not 
input  for  each  problem,  the  previous  title  will  be  used. 

58.  Units  (optional’)  . 

a.  [LN]  •UNITS'  I ! NSTH  TIME 

b.  Definitions: 


’UNIT' 


ran  !  word 


LENGTH  -  ’Vo-character  description  of  length  used 
('FT'  is  tiie  default  value)* 

TIME  -  Throe-character  description  of  time  used 
'  is  the  default  value) 

Note:  un  I  y  r  I  a:  id  *>!'  are  valid  units  of  length.  If  'FT' 

is  used,  then  ' !  i  '  is  the  default  unit  of  force.  If  'M' 
is  used,  then  ’ >  >  '  is  the  default  unit  of  force. 


59.  Water  description. 

a.  [LN]  'WATER'  i 

b.  Definitions: 


MDWT  I!  FLOW 


’WATE'  -  command  word 


coefficient  ot  permeability 

I. eight  of  headwater  against  structure 

change  :r  head  (between  headwater  and  tailwater) 


i  r ec.  t  i  on  (  '  R ' 


to  tiie  right,  'I.'  -■  to  the 


I  C :  ) 


*  A  table  of  factors  for  e. 
(metric)  units  is  present- 


if.  ing  non- SI  units  of  measurement  to  SI 


n  page  J . 


ml  -£>1 


c_.  Note:  Flow  is  a  one-character  description  of  the  direction 
of  flow.  'R'  indicates  flow  to  the  right,  and  'L'  indicates 
flow  to  the  left. 


Fragment  type  descriptions. 


.  'FRAGMENT'  I FACT  DIMENSIONS 


.  Definitions: 

'FRAG*  -  command  word 


IFACT  -  fragment  type 

DIMENSIONS  -  the  dimensions  required  for  the  particular 
fragment  type  chosen  (see  paragraph  61). 

c_.  Note:  A  maximum  of  10  fragments  may  be  entered.  All  frag¬ 
ments  must  be  entered  in  the  direction  of  flow. 


Definitions: 

L  -  The  length  of  the  structure 

A  -  The  height  of  the  confined 
flow  region 


b.  TYPE  2 

'FRAG'  2  T  S 


Definitions: 

T  -  The  height  of  the  flow  region 

S  -  The  distance  that  the  struc¬ 
ture  is  embedded  into  the 
soil 
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£•  TYPE  3. 

'  FRAG’  1  B  T  S 


Def  in  i  t  ions  : 

B  -  The  length  of  the  structure 
resting  on  the  soil 

T  -  The  height  of  the  flow  region 

S  -  The  distance  that  the  struc¬ 
ture  is  embedded  into  the 
soil 


d.  TYPE  4. 
' FRAG ' 


S  FRAGMENT  DIRECTION 


Def init ions : 

B  -  The  length  of  the  structure 
resting  on  the  soil 

T  -  The  height  of  the  flow  region 

S  -  The  distance  that  the 

structure  is  embedded  into 
the  soil 


FRAGMENT  DIRECTION  -  1  or  0 
Indicate  fragment  direction  by  either  a  1  or  a  0. 


Orientation  l 


Orientation  0 


e.  TYPE  5. 


' FRAC’  5  L  T  S 


TvPE  S 


Def initions: 

L  -  The  length  of  the  structure 

T  -  The  height  of  the  flow  region 

S  -  The  distance  that  the  structure 
is  embedded  into  the  soil 


£.  TYPE  6. 

'FRAG'  6  L  T  SI  S2 


Note:  Regardless  of  the  flow  direction,  the  fragment  orienta¬ 

tion  does  not  change;  SI  is  at  the  left  side  and  S2  is  at  the 
right  side. 


Definitions : 

L  -  The  length  of  the  structure 
T  -  The  height  of  the  flow  region 

-  The  distance  that  the  structure 
is  embedded  into  the  soil  on  the 
left 

S0  -  The  distance  that  the  structure 
is  embedded  into  the  soil  on  the 
right 


£.  Restrictions:  S,  SI,  or  S2  should  never  be  greater  than  or 
equal  to  T. 

h.  Note :  Each  fragment  is  aligned  by  using  the  bottom  impervious 
boundary  which  should  always  be  horizontal.  The  dimensions  S, 
SI,  and  S2  may  vary  on  both  sides  of  a  sheetpile,  and  the 
dimension  T  may  vary  from  fragment  to  fragment. 
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62 .  Other  optional  commands. 


a . 

[LN] 

'CONVERSATIONAL' 

provides  conversational  mode 

b. 

[LN] 

' NONCONVERSAT IONAL '  - 

cancels  conversational  mode 

c. 

[LN] 

'MENU' 

draws  menu 

d. 

[LN] 

'GRAPH' 

draws  graph  of  each  fragment  input 

e . 

[LN] 

'MEIN" 

enters  input  through  menu 

jf. 

[LN] 

'COMM' 

provides  list  of  commands 

£• 

[LN] 

'EDIT' 

edits  input  data 

h. 

[  LN  ] 

'EXIT' 

terminates  program  execution 

i . 

[LN] 

’ INFO' 

provides  information  about  program 

Note:  All  of  the  above  commands  may  be  entered  anywhere 

prior  to  the  END  command.  A  more  thorough  explanation  of  these 
optional  commands  is  provided  later  in  this  report. 

63.  Termination. 

a^.  1  END 1  -  Identifies  the  completion  of  data  entry  for  a 

particular  problem. 

Note:  As  previously  mentioned,  a  predefined  data  file  may 

contain  several  problems.  The  END  command  must  be  placed  at 
the  end  of  each  individual  problem.  Once  the  final  problem 
containing  an  END  command  is  processed,  the  user  is  given  the 
option  to  make  additional  runs  or  terminate  the  session. 

A'  ' EXIT1  -  Identifies  the  end  of  a  session. 

Note:  Once  the  'EXIT'  command  is  read,  the  program  is 

immediately  terminated.  Normally,  this  command  would  be 
preceded  by  an  F.ND  command. 


Optional  Methods  of  Input 


Prompt  sequence 

64.  If  the  user  chooses  to  enter  data  from  the  terminal,  then  the 
question 

DO  YOU  WANT  TO  ENTER  DATA  USING  A  PROMPTING  SEQUENCE?  YES  OR  NO. 

is  asked.  If  prompting  is  requested,  then  a  sequence  of  questions  will 

prompt  the  user  for  all  pertinent  data.  An  example  follows: 

ENTER  TITLE  (60  CHARACTERS  MAXIMUM) 

=  EXAMPLE  PROBLEM 

ENTER  UNITS, 

LENGTH  ('FT'  OR  'M')  TIME  (3  CHARACTERS  MAX) 

=  FT  SEC 

ENTER  WATER  DATA, 
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PERMEABILITY  HEADWATER  CHANGE  IN  FLOW  DIRECTION 

(L/T)  HEIGHT  (L)  HEAD  (L)  (LEFT  OR  RIGHT) 

=  .01  20  10  R 

DO  YOU  WANT  TO  USE  A  MENU  AND  CROSSHAIRS  TO  INPUT  FRAGMENT  DATA? 

YES  OR  NO. 

=  Y 

NOTE:  The  sequence  may  go  one  of  two  ways.  If  the  answer  is  'YES', 
then  a  menu  will  be  drawn,  and  fragment  data  may  be  entered  by  selecting 
the  desired  fragment  with  the  cross  hairs.  If  the  answer  is  'NO,'  then 
the  user  will  be  prompted  to  input  the  various  fragment  types  desired. 

All  fragments  must  be  entered  in  the  direction  of  flow.  After  all 
fragment  data  are  input,  the  user  is  given  the  chance  to  edit  the  input. 
The  edit  sequence  is  explained  under  the  EDIT  option.  For  more  informa¬ 
tion  about  the  menu  and  terminal  entry  refer  to  the  optional  command  MEIN 
and  the  discussion  of  the  conversational  mode. 

Nonconversational  mode 

65.  If  the  user  decides  to  enter  data  from  the  terminal  using  the 
nonconversational  mode,  then  an  initial  statement  is  printed  providing 
a  brief  explanation  of  the  input  procedure.  The  program  will  then  prompt 
the  user,  and  data  may  be  entered  by  typing  in  a  command  word  and  any 
subsequent  data.  Initially  the  program  is  in  the  nonconversational  mode. 

An  example  follows: 

DO  YOU  WANT  TO  ENTER  DATA  USING  A  PROMPTING  SEQUENCE?  YES  OR  NO. 

-  N 

ENTER  DATA  USING  COMMAND  WORDS.  TYPE  'END'  TO  COMPLETE  INPUT.  TYPE 
’COMM’  FOR  A  LIST  OF  COMMANDS.  TYPE  'INFO'  FOR  MORE  INFO. 

COMMAND? 

=  NAME  TEST  RUN 

COMMAND? 

=  UNITS  FT  SEC 

COMMAND? 

=  WATE  .01  20  10  R 

COMMAND? 

=  FRAG  2  20  5 

COMMAND? 

=  FRAG  2  20  5 

COMMAND? 

=  END 

*****  INPUT  COMPLETE  ***** 

DO  YOU  WANT  TO  CONTINUE  THE  SOLUTION?  YES  OR  NO. 

=  Y 

The  user  may  continue  with  the  solution  of  the  problem  if  desired. 


Conversational  mode 

66.  When  entering  data  from  the  terminal,  the  user  may  follow  a  conver¬ 
sational  style  of  input  by  entering  the  command  CONV.  This  command  may  be 
entered  any  time  a  command  is  requested.  The  command  words  requiring  addi¬ 
tional  data  may  then  be  typed  in  by  themselves,  and  the  program  will  respond 
with  a  message,  like  those  in  the  prompting  sequence,  which  will  indicate 
the  required  data  needed  for  that  particular  command.  An  example  follows: 

COMMAND? 

=  CONV 

COMMAND? 

=  NAME 

ENTER  TITLE  (60  CHARACTERS  MAXIMUM) 

=TF,ST  RUN 

COMMAND? 

=  WATER 

ENTER  WATER  DATA, 

PERMEABILITY  HEADWATER  CHANGE  IN  FLOW  DIRECTION 

(L/T)  HEIGHT  (L)  HEAD  (L)  (LEFT  OR  RIGHT) 

=  .01  20  10  R 

COMMAND? 

=  FRAGMENT 

ENTER  THE  TYPE  OF  FRAGMENT  I 
=  2 

ENTER:  T  S 
=  20  5 

ENTER  THE  TYPE  OF  FRAGMENT  2 
=  2 

ENTER:  T  S 
=  20  5 


ENTER  THE  TYPE  OF  FRAGMENT  3 

=  (CR)  [If  no  more  fragments  are  to  be  entered,  then  a  carriage  return 
should  be  entered  to  continue.  The  user  may  enter  up  to 
10  fragments. ] 

COMMAND? 

=  END  [Indicates  that  data  entry  is  complete.] 

The  conversational  style  of  input  may  be  switched  off  at  any  time  by  typing 
in  NONC. 


Entry  from  data  file 


67.  Entry  from  a  data  file  is  in  the  same  format  as  the  nonconversa- 


tional  mode  except  that  all  lines  of  input  must  be  preceded  by  a  line 


number : 


30 


10 

NAME 

TEST  RUN 

Enter  title 

20 

UNIT 

FT  MIN 

Sets  units 

30 

WATE 

.01  20  10 

R  -  Enter  water  data 

40 

FRAG 

2  20  5 

Enter  fragment  data 

50 

FRAG 

2  20  5 

60 

END 

This  command  shows  that  input  is 

complete  and  directs  the  program  to  con- 

tinue  with  the  solution. 

Output  Options 

Modification  and  rerun  capability 

68.  After  output  is  complete,  the  question 

DO  YOU  WANT  TO  MODIFY  CURRENT  DATA  AND  RERUN?  YES  OR  NO. 

is  asked.  If  the  user  selects  to  modify  the  data,  then  a  sequence  exactly 
like  the  edit  option  (paragraphs  74-76)  will  be  followed. 

Plotting  of  water  pressures 

69.  If  no  errors  in  the  user's  data  are  found,  the  question, 

DO  YOU  WANT  TO  PLOT  THE  WATER  PRESSURE?  YES  OR  NO. 

is  asked.  If  the  user  responds  with  a  yes,  a  plot  is  drawn  which  includes  a 
diagram  of  the  structure  and  the  resulting  pressure  prisms.  All  pertinent 
dimensions,  water  levels,  and  soil  surfaces  are  shown. 

Use  of  Optional  Commands 


MEIN  command 

70.  When  this  command  is  entered,  the  program  is  automatically  set  for 
the  conversation  mode  as  described  in  paragraph  66.  This  command  must  be 
entered  before  fragment  data.  When  the  command  word  FRAG  is  entered  the 
screen  will  erase,  and  a  menu  of  the  fragment  types  will  be  displayed.  A 
fragment  type  may  then  be  chosen  by  positioning  the  crosshairs  within  the  box 
containing  the  fragment  desired  and  entering  any  character  and  a  carriage 
return.  The  needed  data  for  that  fragment  type  will  then  be  requested. 

After  the  data  are  entered,  the  crosshairs  will  appear  again.  When  fragment 
data  are  complete,  enter  the  character  'E'  to  return  to  the  command  level. 
After  all  input  is  complete,  enter  the  command  END  to  continue  with  the 
analysis.  An  example  is  shown  in  Figure  13. 
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MENU  command 

71.  This  command  will  cause  the  menu  to  be  drawn.  It  may  be  used  to 
help  the  user  distinguish  the  dimensions  of  the  fragments  when  entering  from 
the  terminal,  and  it  may  be  used  later  for  reference  material.  Refer  to 
paragraph  70  for  an  example  of  the  menu. 

GRAPH  command 

72.  This  command  must  be  entered  before  fragment  data.  If  this  command 
is  used  with  the  command  CONV,  a  graph  of  the  fragment  type  entered  will  be 
drawn,  and  the  required  information  will  be  requested.  If  this  command  is  used 
without  the  CONV  command,  then  after  a  fragment  type  and  data  have  been 
entered,  a  graph  of  the  fragment  type  will  be  drawn.  Examples  of  both 
sequences  are  shown  below. 

Sequence  using  CONV  command: 

COMMAND? 

=  CONV 

COMMAND? 

=  GRAPH 

COMMAND? 

=  WATER 

ENTER  WATER  DATA, 

PERMEABILITY  HEADWATER  CHANGE  IN  FLOW  DIRECTION 

(L/T)  HEIGHT  (L)  HEAD  (L)  (LEFT  OR  RIGHT) 

=  .01  20  10  R 

Note:  If  a  prompt  for  only  the  f ragmen*,  data  is  desired,  then  the 

CONV  command  may  be  entered  after  all  other  data  have  been  entered 
but  before  the  FRAGMENT  command. 

COMMAND? 

=  FRAG 

ENTER  THE  TYPE  OF  FRAGMENT  1 
=  •2 

Note:  The  screen  is  erased,  fragment  type  2  is  drawn,  and  the 
required  data  are  requested  (see  Figure  14). 

ENTER  T,  S 
=  20  5 

ENTER  THE  TYPE  OF  FRAGMENT  2 
=  [CARRIAGE  RETURN] 


COMMAND? 


Figure  14.  Computer  response 


Sequence  without  using  CONV  command: 

COMMAND? 

=  GRAPH 

COMMAND? 

=  WATE  .01  20  10  R 

COMMAND? 

=  FRAG  2  20  5 

Note:  The  screen  is  erased  and  fragment  typo  2  is  drawn. 


COMMAND 

73.  After  the  user  enters  this  command,  a  list,  of  the  command  words  will 
appear  at  the  terminal.  It  may  be  entered  any  time  a  command  is  requested. 
Following  is  a  list  of  all  the  commands: 

COMMANDS : 

CONV  -  PROVIDE  CONVERSATIONAL  STYLE  OF  INPUT 
NONC  -  CANCEL  CONVERSATIONAL  STYLE  OF  INPUT 
MENU  -  DRAW  MENU 

GRAP  -  DRAW  GRAPH  OF  FRAGMENT  TYPE  AND  INPUT  DATA 
MEIN  -  DRAW  MENU  AND  INPUT  FROM  CROSSHAIRS 


FRAG  -  ENTER  FRAGMENT  TYPE  AND  DATA 

WATE  -  ENTER  WATER  DATA 

END  -  INPUT  COMPLETE 

COMM  -  PRINT  LIST  OF  COMMANDS 

EDIT  -  EDIT  INPUT 

EXIT  -  EXIT  PROGRAM 

UNIT  -  ENTER  LENGTH  AND  TIME  UNITS 

INFO  -  PROVIDE  INFORMATION  ABOUT  PROGRAM 

EDIT  command 

74.  Entering  this  command  invokes  a  sequence  whereby  the  input  data 
may  be  changed.  It  is  the  same  sequence  used  in  the  Prompting  Sequence  and 
the  Modification  of  Data  Options.  If  data  are  being  entered  from  the  terminal, 
the  EDIT  command  may  be  given  once  the  input  is  finished.  Input  is  assumed 

to  be  complete  upon  completion  of  the  edit  sequence. 

75.  The  edit  sequence  is  in  four  sections:  the  title,  the  units,  the 
water  data,  and  the  fragment  data.  Any  or  all  of  these  sections  may  be 
edited.  In  addition,  the  fragment  section  has  several  commands  to  change 
fragment  data.  They  are  as  follows: 

a.  CHANGE  -  This  command  is  used  to  change  fragment  types. 

The  new  data  needed  will  be  requested. 

b.  DELETE  -  This  command  will  delete  a  specified  fragment. 

£.  ADD  -  This  command  will  add  a  fragment  to  the  existing 

arrangement  of  fragments  at  the  location  specified. 

The  CHANGE,  DELETE,  and  ADD  commands  are  set  up  in  the  following  manner: 

OPTION  A  B 

CHANGE  FRAGMENT  NO.  FRAGMENT  TYPE 

DELETE  FRAGMENT  NO. 

ADD  FRAGMENT  NO.  FRAGMENT  TYPE 

The  user  should  enter 

OPTION  A  B 

For  example, 

CHANGE  3  2 

DELETE  4 

76.  A  detailed  description  of  each  command  follows: 

a.  CHANGE  command. 

OPTION  =  CHANGE 

A  =  Number  of  the  fragment  to  be  changed. 

B  =  New  fragment  type  desired. 

After  this  command  is  entered  the  information  required  tor 
the  new  fragment  type  will  be  requested. 


4,  5,  and  6  can  only  be  located  within  the  structure;  they  cannot  act  as 
entry  or  exit  fragments.  Finally,  a  maximum  of  10  fragments  may  be  entered 
for  each  problem. 

Notes  on  error  messages 

81.  In  addition  to  the  error  messages  previously  discussed,  other 
error  messages  may  be  encountered  in  the  execution  of  the  program.  Two 
example  conditions  which  will  cause  error  messages  are:  (a)  insufficient 
number  of  values  entered  for  the  specified  command  word  and  (b)  an  incorrect 
parameter  detected  in  the  data  list. 

82.  To  correct  either  of  the  above  errors  while  using  the  prompt 
sequence  or  the  conversational  mode,  the  user  is  required  only  to  re-enter 
the  variable  list.  The  user  will  not  be  allowed  to  continue  until  the 


required  data  are  entered. 

83.  In  the  nonconversational  mode  or  while  entering  data  from  a  data 
file,  the  user  must  re-enter  the  command  word  along  with  the  required  data. 
The  user  may  select  to  re-enter  the  command  in  which  the  error  was  detected 
or  continue  with  another  command.  The  only  exception  to  this  rule  occurs 
when  data  are  specifically  asked  for  again. 

84.  In  addition,  if  entry  is  from  a  file  and  an  error  message  is 
encountered,  several  options  may  be  exercised:  (a)  the  particular  command 
may  be  corrected,  (b)  another  command  may  be  entered,  or  (c)  a  carriage 
return  may  be  entered  which  will  result  in  the  next  line  of  data  being  read 
from  the  file  and  the  incorrect  line  being  ignored. 


APPENDIX  A:  EXAMPLE  COMPUTER  RUNS  AND  HAND 
VERIFICATIONS 


1.  This  appendix  contains  example  runs  and  hand  computations.  Solutions 
are  provided  for  typical  problems  that  can  be  analyzed  by  the  method  of 
fragments.  The  examples  are  intended  to  show  how  the  various  fragment  types 
can  be  utilized  to  model  a  flow  region  and  how  the  form  factor  for  each 
fragment  type  is  calculated.  The  hand  computations  are  provided  to  verify 
the  results  from  CFRAG. 

2.  In  the  diagrams  accompanying  each  problem,  the  different  fragments 
are  identified  by  large  circled  numbers,  and  points  of  pressure  calculations 
are  identified  by  small  circled  numbers. 

3.  To  analyze  this  structure,  the  user  must  divide  the  flow  region 
shown  into  three  fragments.  The  analysis  is  performed  from  left  to  right 
corresponding  to  the  direction  of  flow. 


Example  Al:  Flow  Under  Simple  Structure  Without  Sheetpiles 


Fragment  number  1 

4.  Fragment  number  1  is  an  entrance  condition  and  is  a  type  2  fragment. 


S  = 

T  = 

<J>  = 

m  =  sin  2(2  0 )  =  °-3090 

m2  =  0.0955 

(m')2  =  1  -  m2  =  0.9045 
From  tabulated  values  of  elliptic  integrals  of  the  first  kind, 

K  =  1.610 
K'  =  2.600 
|r  =  0.619 

$  =  0.619 
Fragment  number  2 

5.  Fragment  number  2  is  a  type  1  fragment. 

L  =  18 
A  =  16 


20 

K 

K' 


,  .  irS 

where  m  =  sin  — 


18 

16 


=  1.125 

Fragment  number  3 

6.  Fragment  number  3  is  a  type  2  fragment  and  is  an  exit  condition. 


S  =  4 
T  =  2<- 


$ 


K  ,  ,  ttS 

where  m  =  sin 


A2 


7.  This  fragment  is  identical  with  fragment  number  1;  therefore,  the 
form  factor  is  also  the  same. 


*  =  0.619  . 

Head  loss  calculations 

8.  Head  loss  calculations  are  as  follows: 
£*  =  2(0.619)  +  1.125  =  2.363 


h 

m 


h4>m 


14(0.619) 

2.363 


3.67 


14(1.125) 

2.363 


6.67 


h3  =  3.67 


Total  head  loss  =  3.67  +  6.67  +  3.67  =  14.0  ft 
Pressure  calculations 

9.  Pressure  loss  calculations  are  as  follows: 


z„ 

-  h. 

2 

L 

. 

=  16(62.4)  =  998.4  psf 
P2  =  (20  -  3.67)  62.4  =  1019.0  psf 
P3  =  (20  -  (3.67  +  6.67))  62.4  =  603  psf 
P4  =  (16  -  14)  62.4  =  124.8  psf 

Total  flow  calculation 


9  =  2  363  ~  =  85.32  ft^/day  per  foot  of  width 

Exit  gradient  calculation 


h3ir 

2KTm 


_ 3.67(tt) _ 

2(1. 610) (20) (0. 3090) 


0.579 


A3 


CONFINED  FLOW  -  METHOD  OF  FRAGMENTS 
TINE  IS  3S  5S 


DATE  •  8/  8 ■'83 


TITLE  -  Sliplt  Strvctor*  mtkttl  SknipllM 
0  •  8S  2829  (FTM2/DAV) 

K  •  14  4808  (FT'DAV ) 

O'K  •  S  92  (FT ) 

TOTAL  HEAD  LOSS  •  14  ••  (FT) 


FRAG 

NO 

FRAC 

TYPE 

L 

(FT) 

A 

(FT) 

B 

(FT) 

T 

(FT) 

SI 

(FT) 

S2 

(FT) 

FORM 

FACTOR 

1 

2 

28  88 

4  88 

8  62 

2 

1 

18  88 

16  88 

1  13 

3 

2 

28.88 

4.88 

8.62 

EXIT  GRADIENT  •  •  5798 


I 

RESULTANT  FORCES  ON  STRUCTURE 


LATERAL  FORCE  UPLIFT  FORCE 
HEADUATER  SIDE 

(LBS)  (LBS) 


LATERAL  FORCE 
TAILUATER  SIDE 
(LBS) 


12822.2 


14681  6 


DO  YOU  WANT  TO  PLOT  UATER  PRESSURES?  YES  OR  NO 


1581  8 


Figure  Al.  Program  output  for  Example  A1 


IN  POUNDS 


10.  To  analyze  this  structure,  the  user  must  divide  the  flow  region 
shown  into  two  fragments.  The  analysis  is  performed  proceeding  from  left 
to  right. 


Fragment  number  1 

11.  Fragment  number  1  is  a  type  2  fragment  and  an  entrance  condition 


m  =  sm 
2=  0.50 


'  TT  (  10) 
2(20) 


=  0.7071 


m 


m'2  =  1-m2  =  0.50 


S  =  10 
T  =  20 

$  =  — 7-  where  m  =  sin 
K 


ttS_ 

2T 


From  tabulated  values  of  elliptic  integrals  of  the  first  kind. 


Fragment  number  2 

12.  Fragment  number  2  is  a  type  3  fragment  and  an  exit  condition. 


b  =  16 
T  =  17 


K  ,  ttS  */  ,2  irb  2  tt S 

•jrr  where  m  =  cos  —  tan  h  —  +  tan  — 


m  =  cos 


m  =  0.9383 


m“  =  0.8803 


tan  h 


/tt(16)  \  2  /  tt (7)  \ 

y (2(17)  )  tan  y  2(17) J 


From  tabulated  values  of  elliptic  integrals  of  the  first  kind, 


£r  =  1.54 


4>  =  1.54 


Head  loss  calculations 


13.  Head  loss  calculations  are  as  follows: 


=  1.00  +  1.54  =  2.54 


h  =  h^m 

m  IT 


h  =  10(1.00)  =  3.94  ft 
2.54 


h2  =  101^54!  __  6>06  ft 
2.54 

Head  loss  gradient  in  fragment  2, 

Gradient,,  =  -J>~ =  0.2635 

Total  head  loss  =  3.94  +  6.06  =  10.0  ft 
Pressure  calculations 

14.  Pressure  loss  calculations  are  as  follows: 

p“N  _hJ  y“ 

?1  =  15(62.4)  =  936  psf 

?2  =  (25  -  3.94)  62.4  =  1314.1  psf 


A7 


P3  =  [18  -  (3.94  +  0.2635(7))!  62.4  =  762.5  psf 
P4  =  (18  -  (3.94  +  6.06))  62.4  =  499.2  psf 

Total  flow  calculation 


CONFINED  FLOW  -  METHOD  OF  FRAGMENTS 

Tint  16  38  42  DATE  8/  8/83 

TITLE  -  Ue  L  r  with  Upsc'eaa  Shaatplla 
0  •  56  7307  (FTM2/DAV) 

K  •  14  4000  (FT/DAY > 

0/K  •  3  94  (FT) 

TOTAL  HEAD  LOSS  •  10  00  (FT) 


FRAG 

FRAG 

FORM 

HEAD 

NO 

TYPE 

t 

A 

B 

T 

SI 

S2 

FACTOR 

LOSS 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

1 

2 

20  00 

10  00 

1  00 

3  94 

2 

3 

16  00 

17  00 

7  00 

1  54 

6  06 

WARNING  -  THERE  IS  NO  EMBEDMENT  ON  THE  TAILUATER  SIDE 
EXIT  GRADIENT  IS  INFINITE  AND  PIPING  NAV  OCCUR 


RESULTANT  FORCES  ON  STRUCTURE 

LATERAL  FORCE  UPLIFT  FORCE  LATERAL  FORCE 

HEADWATER  SIDE  TAILUATER  SIDE 

(LBS)  (LBS)  (LBS) 


18270  8  10091  8  9264  3 


DO  YOU  WANT  TO  PLOT  UATER  PRESSURES?  YES  OR  NO 
-Y 


Figure  A3.  Program  output  for  Example  A2 


A  8 


0.879 


F  ragi nor.  t  ntiRihe  r  - 1 

19.  Fragment  niu::!-**  i  *  .  •  an  ex  i  t  rendition  .aid  is  type  fragment 


( ■ (14) 

V  ""  V’TJ,7V 

m~  =  11.9603 


0.  .69 


n..; 


wno  re  m  -  sm  -- 


From  tabulated  values  of  elliptic  integrals  of  the  tirst  kind, 
K  =  1.908 
K  ’  =  1.806 

fr*  l'0h 

1  =  1.06 


Head 

1  os- 

cal  c u 

1  a  t  i ons 

ll). 

Head 

loss  calculations  are 

as  follow.* 

;  ;  = 

0.9! 

+  0. 

391  4  0.8/9  +  1.06  -  : 

3.43.’ 

1,  = 
n 

him 

hl  ’ 

i  9  (  . 
3. 

,  'h:,.  1 

.  'I  ' 

7  .  h 

V  = 

3  9  f  . 

3 . 

■)“]) 
'Vi  7  ' 

(•rad  ient  .. 

■  v  Wt  “  "•ml 

"i  = 

1  9  ( . 
~  "3". 

,  8  7  9 ) 

,9  3..’ 

4 . 8  7 

Orad  ient 

.3 

-  vvidVis  =  °-n91 

19(1.  Ur,  , 


Total  head  loss  *  .  * .  .H  i  *  9. 67  -  1‘i.ti  ft 


Pressure  calculat ions 


21.  Pressure  loss  r aleul ations  are  as  follows: 


- 

“ 

p  = 

Z2 

-  hL 

. 

P  =  24(62.4)  1497.6  psf 

P,,  =  (36  -  5.0)  62.4  =  1934.4  psf 

P3  =  [28  -  (5.0  +  0.1211(8.0))]  62.4  =  1374.7  psf 

P4  =  [28  -  (5  +  0.1211(19))]  62.4  =  1291.6  psf 

P5  =  [36  -  (5  +  3.27)]  62.4  =  1730.4  psf 

P£  =  [28  -  (5  +  3.37  +  .1391(8))]  62.4  =  1161.7  psf 

o 

?7  =  [28  -  (5  +  3.27  +  .1391(23))]  62.4  =  1031.5  psf 

Pg  =  [40  -  (5  +  3.27  +  4.87)]  62.4  =  1676.0  psf 

P9  =  7(62.4)  =  436.8  psf 
Total  flow  calculation 

kh 

q  1$ 

q  =  =  79.72  ft^/day  per  foot  of  width 

Exit  gradient  calculation 

=  h47  _  _ 5.87(0 _ 

iE  2KTm  2(1 .908) (26)  (0.7485) 


0.248 


CONFINED  FLOW  -  NETHOO  CF  FRAGMENTS 
TIME  16  43  4X 


DATE  8/  8 ''83 


TITLE  -  Daw  with  Hultlplo  Shaatpl <as 

0  ■  79  8086  (FTIia/DAV) 

K  -  1 4  4000  iFT/DAV) 

O/F  •  S  E 4  (FT) 


TOTAL  HEAD  LOSS  ■  19  @0  <  F  T ) 


FRAG 

MO 

FRAG 

TYPE 

l  A 

<FT)  <fr> 

8 

( F  T  ) 

T 

(FT  ) 

SI 

(FT  ) 

sa 

(FT) 

FORT! 

FACTOR 

HEAD 

LOSS 

(FT) 

1 

a 

38  00 

13  00 

0  90 

S  00 

a 

s 

1  1  00 

34  00 

8  00 

0  59 

3  87 

3 

6 

IS  00 

34  00 

8  00 

ia  00 

0  88 

4  87 

4 

a 

36  00 

14  00 

1  06 

S  B6 

EXIT 

GRADIENT 

•  0  3477 

RESULTANT  FORCES  CM  STRUCTURE 


LATERAL  FORCE  UPLIFT  FORCE 

HEADUATER  SIDE 

<LBS1  (LBS) 


LATERAL  FORCE 
TAILUATER  side 
<  LBS  ) 


J8563  5 


31109  1 


16317  5 


DO  VOU  UANT  TO  PLOT  WATER  F>PESSUR£S?  YES  OR  HO 

-V 


•  i;  >  ,»n,  i  i :  ‘  :  i :  *  fur  Ky.'implu  A3 


D*«  Ml  ill  Hw  ttijpU  SKwtBllM 


Figure  A6.  Plot  of  pressures  for  Example  A3 


i.xample  A3 .  Structure  with  a  Sheetpile  and  a 
"Step-down"  in  the  Base 


22.  To  analyze  this  structure,  the  user  must  divide  the  tlow  region 
shown  into  four  fragments.  The  analysis  is  performed  from  left  to  right. 


i 

tx: 


©  i 


t  r.ifcr.ont  mi:::: 


23.  Frarrru  r.  r. 


.a,  entrance  condi t  ion  and  is  a  type  2 


t  ragr.ent . 


where  m  =  sin 


Fragment  number  4 

26.  Fragment  number  4  is  an  exit  condition  and  a  type  2  fragment 


S  =  10 
T  =  35 


*  K  .  ttS 

•P  =  ^rr  where  m  =  sm 


V sln  (fm)  *  °-4339 

m“  =  0.1883 

m'2  =  1  -  m2  =  0.8117 


From  tabulated  values  of  elliptic  integrals  of  the  first  kind. 


K  = 

K*  = 
K_ 

K' 


1.654 

2.285 

0.724 


t  =  0.  724 


Head  loss  calculations 

27.  Head  loss  calculations  are  as  follows: 


T.<t  =  0.546  +  0.405  +  0.823  +  0.724  =  2.498 
,  him 

h  = 

m 


-  20(0.546) 

1  2.498 
=  20(0.405) 

2  2.498 
_  20(0.823) 

"3  2.498 

=  20(0. 724) 

4  2.498 

Total  head  loss 


4.37 

3.24 

6.59 

5.80 

=  4.37  +  3.24 


Gradient^ 

Gradient^ 

+  6.59  +  5.80 


3.24 
10  +  10 
6.59 
5  +  20  ' 


0.1620 

0.2636 


20.0  ft 


Pressure  calculations 

28.  Pressure  loss  calculations  are  as  follows: 

r  —y 


Pj  =22(62.4)  =  1372.8  psf 

P2  =  (27  -  4.37)  62.4  =  1412.1  psf 

P3  =  [27  -  (4.37  +  0.1620(10))]  62.4  =  1311.0  psf 

P4  =  (37  -  (4.37  +  3.24))  62.4  =  1833.9  psf 


A18 


P5  =  (32  -  (4.37  +  3.24  +  .2636(5))]  62.4  =  1439.7  psf 
P6  =  [32  -  (4.37  +  3.24  +  6.59)]  62.4  =  1110.7  psf 

P7  =  2(62.4)  =  124.8  psf 

Total  flow  calculation 


kh 

q  =  1$ 


14.4(20.0)  llr  „„  ^ . 3 ,  ,  ,  .  .... 

q  =  -  y  -  =  115.29  ft  /day  per  foot  of  width 


Exit  gradient  calculation 


_  h4Tr  _  _ 5.80(tt) 

E  2KTm  2(1. 654) (35) (0. 4339) 


0.363 


CONFINED  FLOW  -  NETMOD  OF  FRAGNENTS 

TINE'  164148  DATE'  8/8/83 


TITLE  -  Dm  «lth  Interior  Sheetpile  with  Non-Uni  fore  Bene 

0  •  IIS  2679  (FTM2/DAY) 

t  •  14  4888  (FT/DAV) 

O/tC  •  8  88  (FT) 

TOTAL  HEAD  LOSS  •  28  88  (FT) 

FRAG  FRAG 
NO  TYPE 


1 

2 

3 

4 


L 

A 

B 

T 

SI 

S2 

FACTOR 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

3S  88 

5  88 

8  55 

18  88 

38  88 

18  88 

8  41 

28  88 

25  88 

5  88 

8  82 

35  88 

18  88 

8.72 

EXIT  GRADIENT  •  8  3624 


RESULTANT  FORCES  ON  STRUCTURE 


LATERAL  FORCE  UPLIFT  FORCE 
HEADWATER  SIDE 

<LBS)  (LBS) 


LATERAL  FORCE 
TAILUATER  SIDE 
(LBS) 


37784  2  39118.1 


14482  4 


DO  YOU  WANT  TO  PLOT  UATER  PRESSURES?  YES  OR  NO 


Figure  A7.  Program  output  for  Example  A4 


Paa  Mtih  Intarlor  SHaat 


Figure  A8.  Plot  of  pressures  for  Example  A4 


Ai'i'KXiu  \  :  •  -  'Mi’AK  i  SON  ANALYSIS  MKT  HODS 

1.  In  '.his  appendix,  exampl os  are  presented  to  compare  the  method  of 
fragment  w  •»  ot  h-->-  -onraon  ’  v  "s  •«  net  hoc's  ->f  ana  !  vz  ing  se  ■png*.  Included  in 
the  >■!  '!”p.i  ■'  i  sot;  a  re  so.iit  ions  ci-mg  .re  me  thou  of  creep,  a  l  ini  to  element 
met  hee  ,  and  ' outlets.  The  f  in  it  e  element  program  is  titled  SEEPKG,  and  it 
was  written  hv  Mr.  Fred  Id"  -v,  ADP  Center,  Waterways  Experiment  Station.  The 
ye  :::  i.  u  •  ,ie  COIL'S  )"a  v  ,  XSl'C l. 

i .  T  ■•*.>'  cnl :  ft  force  c  mi  tho  location  of  the  resultant  as  well  as 
i  ;,i  total  t  to.v-  n:u!  '<it  grad  iort  a*u  -an  pa"  -c!  lor  each  problem.  A  plot  of 
tiie  pressures  obt  \  '•  -  vd  hv  c  ci  r  i  n  d  0,  included  as  a  visual  :  id  in  the 
compel'  ison . 


Exanp  *  t  A  Single^  Sheetpile 


3  .  The  p  rob  1  ■  n  '  !  1  r.  ! 
a  single  sheetpile. 


’  .re  31  is  to  analyze  seep.. go  around 


•i 


4.  The  given  problem,  modeled  by  two  type  2  i  r.  s  re?'  t ,  m.r  be  solved 

by  the  method  of  fragments  as  shown  in  Figure  Hie.  The  re-suit:.  from  CFRAC  are 
shown  in  Figures  Bib  and  Blc. 

5.  The  solution  of  the  problem  by  a  ilownet  is  shown  in  Figure  Bid. 

The  grid  used  in  the  finite  element  analysis  is  shown  in  Figure  Ble.  The 
sheetpile  modeled  by  a  1-ft-wide  gap  is  shown  in  Figure  Ble. 

6.  Pressure  diagrams  obtained  by  each  method  are  compared  in  Figure  Blf. 
Uplift  and  total  flow  and  exit  gradients  determined  by  the  three  methods  are 
compared  in  Tables  Bla  and  Bib,  respectively. 


B2 


Shaatpl !• 


Figure  Blc.  Plot  of  pressures  for  Example  Bl,  method  of  fragments 


Figure  Bid.  Flownet  analysis  of  single  sheetpile 


-0.1320E  02 
0.13-42E  03 
-0.1220E  02 
0.7320E  02 


Figure  Ble.  Grid  for  finite  element  of  single  sheetpile 


Uplift  Comparison,  Single  Sheetpile 


Lateral  Force  on  Location  of  Re-  %  Difference  from 

Headwater  Side  sultant  from  Top _ Finite  Element 


Method 

lb 

ft 

Lateral  Force 

Resultant  ' 

MOF 

50668.8 

27.512 

2.55 

0.670 

CREEP 

50668.8 

27.512 

2.55 

0.  670 

FLOWNET 

53272.0 

27.924 

2.46 

0.816 

FEM 

51979.2 

27.697 

Table  Bib 

Total  Flow  and  Exit  Gradient  Comparison, 
Single  Sheetpile 


%  Difference  from 

Total  Flow  Finite  Element 


Method 

f  t-Vday 

Exit  Gradient 

Flow 

Exit  Gradient 

MOF 

244.8 

0.4338 

9.16 

11.29 

FLOWNET 

240.0 

0.37 

11.14 

27.04 

Example  B2  :  Dam  with  Sheetpile  at  Toe 

7.  The  problem  illustrated  in  Figure  B2  is  a  typical  example  requiring 
sound  engineering  judgment  to  properly  model  a  flow  region. 


8.  The  problem  was  divided  into  three  fragments  (Figure  B2a)  by  assum¬ 
ing  that  the  upstream  and  downstream  keys  have  negligible  thicknesses.  These 
simplifications  were  necessary  in  order  to  obtain  rational  results.  From  an 
earlier  discussion  it  was  determined  that  the  basic  assumption  in  the  develop¬ 
ment  of  the  method  of  1 ragments  was  that  the  equipotential  lines  must  be 
approximated  by  vert  teal  1 itu-s.  Had  the  geometry  of  the  structure  not  been 
simplified,  a  flow  region  with  six  t ragments  (Figure  B2b)  would  have 
resulted.  Notice  fii.it  t  hose  equipotential  lines  do  not  approximate  the 
actual  equipotential  1  in. s  shown  in  Figure  B2a.  Consequently,  results  using 
this  division  would  he  inon  rt'ct  . 

9.  The  results  ot  tin  program  "un  are  shown  in  Figures  B2c  and  B2d , 
and  the  flownet  for  this  problem  is  shown  in  Figure  I12e.  The  grid  used  in 


the  finite  element  analysis  is  shown  in  Figure  B2f.  The  plot  of  the  pres¬ 
sure  prism  obtained  by  each  method  is  shown  in  Figure  B2g.  Results  obtained 
using  each  method  are  compared  in  Tables  B2a  and  B2b. 
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C ONF I ME D  FLO U  -  METHOD  Of  FRAGMENTS 
rrne  16  S3  46 


DATE  8*  8/83 


TITLE  -  Da*  with  Sheets*  1»»  *  L  Toe 
Q  •  ?2  SS?3  (FTII2' Dav  i 

K  -  14  4*99  (FT.'DAV) 

Q/K  •  5  04  (FT) 

TOTAL  HEAD  LOSS  -  IS  09  <f  Ti 


FRAG 

MO 

frag 

TVPE 

l 

.  f  T  ■ 

A 

(  KT  ) 

8 

1 E  T  ) 

T 

'FT) 

El 
t  -  T 

SP 
!  f  T  1 

P0Rfl 
f AC  TOR 

HE* 

L01 

(PI 

1 

?. 

1 4  ee 

13 

0  5« 

2 

2 

6 

Si  C-0 

63  S3 

f.  00 

32  S3 

1 

7 

3 

P 

ee 

39  ?0 

1  32 

5 

EXIT  GRADIENT  •  P  <?*?’ 


Pf SULTAN?  FOR  f  m  STRUv;TURL 


»  Avi  »AL  r  ORCf  UP!.  IK'  '  ORCfc 

H'AruATKR  OI14: 

iLBS)  <LHS) 


L  A  ’  *  .V  A  »  *  "  b 

TATL^ATfP  - . 

;  LB  3 


28  n:  4  < 


10W3  9 


61779  8 


DO  VOU  UA«T  rC  P !'.«.  T  .,  rV  *■»',  -iXfr ',v  vf  0  -)R  ho 


xhin  • 

XPAX  ■ 
xniN  • 
vpirtx  • 


10  FT.-  f. 


'A/W 


► 

r 

M; 


i 

i 

-j- 


~7f~T- 


“7n 


.  V ;l.. 


/'■X/Wt 


\y 


\/\y- 


Figure  B2e. 


F 'li  win  t  analysis  of  Kxamp  1  e  B2 :  dam 
v:tn  a!: net  pi  j  e  at  t<'t 


-0.2008E  03 
0.2888E  03 
-0.4C80E  02 
0.1143E  03 


Figure  B2f.  Orid  u.-.<-d 


:  or 


finite-  element  analysis  of  Example  B2 


Total  Uplift  Location  oi  kesultant  _ Finite  Element 


Method 

Force,  lb 

from  i.eft,  ft 

t ’p 1 i ft  Force 

Resultant  Location 

MOF 

106053.8 

40.754 

5.07 

1.73 

CREEP 

115041.6 

40.791 

3.07 

1.64 

FLOWN ET 

112973.3 

4 1 . 589 

1  .25 

0.294 

FEM 

111568.7 

'. !  .  49/ 

Tab Le  Bib 

Total  Flow  and  Fxit  Gradient  Comparison ,  Dam  with 
Sheetpile  at  Toe 


%  Difference 

from 

Totai 

L  Flow 

Finite 

Element 

Method 

ff 

/day 

Exit  Gradient 

Flow 

Exit 

Gradient 

MOF 

72. 

.6 

0.080  7 

0. 

,552 

16 

.70 

FL0WNET 

72. 

.0 

0.1 

0. 

,277 

4 

.71 

10.  Example  P2  was  reworked  with  the  flow  revet-. •■! 


In  essence,  this 


places  the  sheetpile  at  the  heel  of  the  dan  instead  ot  at  the  toe.  the  pro¬ 
gram  output  is  shown  in  Figures  B2h  and  iiJ i .  The  flownet  shown  in  Figure  B2j 
is  merely  the  reverse  of  that  shown  in  Figure  B2f.  The  finite  element  grid  is 
also  the  same  as  shown  in  B2e.  A  plot  ot  the  pressure  prism  tor  each  method 
is  presented  in  Figure  B2k.  The  methods  are  compared  in  Tables  B2c  and  B2d. 

CONFINED  FLOW  -  NETHOO  OF  FRAGMENTS 

TINE  17  3  39  DATE  8/  8/83 


TITLE  -  Daa  with  Shaslpi Is  st  Nest 
0  •  72  5578  (FTM2/DAV) 

K  •  14  4000  (FT/DAV) 

0/K  •  5  04  <FT ) 

TOTAL  HEAD  LOSS  ■  15  00  (FT) 


FRAG 

NO 

FRAG 

TYPE 

L 

(FT) 

A 

(FT  ) 

B 

(FT) 

T 

(FT) 

SI 

(FT  ) 

S2 
(FT  ) 

FORN 

FACTOR 

HEAD 

LOSS 

(FT) 

1 

2 

74  00 

38  00 

1  02 

5  1- 

2 

6 

85  00 

68  S0 

32  50 

5  00 

1  41 

7. 11 

3 

2 

74  00 

10  50 

0  54 

2  7! 

EXIT  GRADIENT  •  0  1657 


RESULTANT  FORCES  ON  STRUCTURE 


LATERAL  FORCE  UPLIFT  FORCE 

HEADUATER  SIDE 

(LBS)  (LBS) 


LATERAL  FORCE 
TAILUATER  SIDE 
(LBS) 


98865  8 


84890  2 


8395  1 


DO  YOU  UAKT  TO  PLOT  UATER  PRESSURES?  VES  OR  NO 
*Y 


Figure  B2h.  Program  output  for  Example  B2,  dam 
sheetpile  at  heel 


V  *  Lwt  OF 


Total  Uplift 

Location  of  Resultant 

Finite  Element 

Method 

Force,  lb 

from  Left,  ft 

Uplift  Force  Resultant  Location 

MOF 

84864.0 

40.313 

10.46 

2.30 

CREEP 

90753.6 

40.189 

3.75 

2.60 

FLOW 

NET 

92712.7 

41.268 

1.62 

0.046] 

FEM 

94226.5 

41.249 

Table  B2d 

Total  Flow  and  Exit  Gradient  Comparison,  Dam  with 
Sheetpile  at  Heel 


Method 

MOF 

FLOW 

NFT 

FFM 


Total  Flow 
f  t^/day 

%  Difference  from 
Finite  Element 

Exit  Gradient 

Flow 

Exit  Gradient 

72.6 

0. 16S7 

0.552 

7.13 

7  2.0 

0.  !  5 

0.277 

2.83 

72.2 

0. J 543 

APPENDIX  C:  PARAMETRIC  STUDY 


1.  A  brief  parametric  studv  was  performed  to  compare  the  method  of 
fragments  witli  the  method  of  creep.  The  accuracy  ol  each  method  was 
evaluated  using  the  finite  element  method  with  the  grid  shown  in  Figure  Cl. 

In  the  analysis,  Location,  length,  and  number  of  sheetpiles  were  varied  while 
the  width  and  depth  of  the  flow  region  were  held  constant. 

2.  The  total  uplift  pressure,  the  location  of  the  resultant  uplift,  the 
total  flow,  and  the  exit  gradient  were  compared.  Also,  a  plot  of  the  uplift 
pressure  diagrams  is  included. 

3.  The  results  show  that  the  assumption  of  vertical  equipotent ial 
lines  between  fragments  made  in  the  method  of  fragments  is  more  accurate 
when  the  length  of  the  sheetpiles  increases  and/or  when  two  sheetpiles  are 
present  under  the  structure.  Corresponding  to  a  decrease  in  the  difference 
between  the  actual  and  assumed  equipotent ial  lines  was  an  increase  in  the 
accuracy  of  the  method  of  fragments  solutions. 


xniN  •  -0.8000E  01 
XHAX  •  0.8800E  03 
YOIN  •  -0.8C00E  01 
vrw  •  0.380OE  02 


Figure  Cl.  Grid  used  in  finite  element  analysis  of 
the  various  dams  in  the  parametric  study 


Cl 


EXAMPLE  6 


EXAM1 


Figure  C2.  (Concluded) 


CONFINED  FLOU  -  METHOD  OF  FRAGMENTS 

TIME  I?  711  DATE-  8x8'83 

TITLE  -  Daa  with  no  Shaetpllas 
0  •  94  2068  (FTM2/DAV) 

K  •  H  4000  (FT/DAV) 

Q/K  •  6  54  (FT) 

TOTAL  HEAD  LOSS  •  10  00  (FT) 


FRAG 

FRAG 

FORn 

HEAD 

NO 

TYPE 

L 

A 

B 

T 

SI 

S2 

FACTOR 

LOSS 

(FT) 

(FT  ) 

(FT! 

(FT) 

(FT) 

(FT) 

(FT) 

1 

3 

10  00 

30  00 

0 

0  76 

5  0 

2 

3 

10  00 

30  00 

0 

0  76 

5  0 

WARNING  -  THERE  IS  NO  EMBEDMENT  ON  THE  TAILUATER  SIDE 
EXIT  GRADIENT  IS  INFINITE  AND  PIPING  MAV  OCCUR 


RESULTANT  FORCES  ON  STRUCTURE 

LATERAL  FORCE  UPLIFT  FORCE  LATERAL  FORCE 

HEADWATER  SIDE  TAILUATER  SIDE 

(LBS)  (LBS)  (LBS) 


3120  0  6240  0  0 


DO  VOU  UANT  TO  PLOT  UATER  PRESSURES?  VES  OR  MO 

-V 


Figure  C3a.  Program  output  for  Example  1 


C3 


n«» 


Figure  C3c.  Plot  of  uplift  pressures  for  each 
method  for  Fxample  1 


CS 


Table  ("3a 

I  plift  Comparison  for  Example  1 ,  Dam  With  No  Sheetpiles 


Total  Dplitt  Location  of  Resultant 


Z  Difference  from 
Finite  Element 


Method 

Force,  lbs 

from  Left  ft 

Uplift  IN 

MOF 

t,J  .  0 

6.667 

0.029 

CREED 

6240.0 

6 . 6b  7 

0. 026 

FEM 

6341  .8 

7.186 

Table  C3b 

fotal  Flow  and  Exit  Gradient  Comparison  For  Example  1, 


Dam  With  No  Sheetpiles 


7,  Difference  from  0M 

. 

Total  FI ow 

Finite  Element  1 

• ,  ~ 

Method 

(ft3 /day) 

Exit  Gradient 

Flow  Exit  Gradient  -  4 

MOF 

94.2 

* 

5.97 

C 

FEM 

100 

* 

5.  As  shown  in  Figure  C2b,  Example  2  is  a  dam  with  one  6-ft  sheetpile 
at  toe.  The  method  of  fragments  solution  is  shown  in  Figures  C4a  and  C4b. 

A  plot  of  the  uplift  pressures  for  each  method  is  shown  in  Figure  C4c.  The 
uplift  comparison  for  Example  2  is  shown  in  Table  C4a,  and  the  total  flow  and 
exit  gradient  are  compared  in  Table  C4b. 


CONFINED  FLOU  -  HETHOD  OF  FRAGMENTS 

TINE  17  14  59  DATE  8/  8/83 

TITLE  -  Dh  uith  ona  6'  ShMtplIl 
0  ■  81  S057  (FTM2/DAY) 

K  •  14  4000  (FT/DAV) 

0/K  •  S  66  (FT) 

TOTAL  HEAD  LOSS  •  10  00  (FT) 


FRAG 

FRAG 

FORfl 

HEAD 

NO 

TYPE 

L 

A 

B 

T 

SI 

S2 

FACTOR 

LOSS 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

(FT) 

1 

3 

20  00 

30  00 

6  00 

1  15 

6  49 

2 

2 

30  00 

6  00 

0  62 

3  SI 

EXIT  GRADIENT  •  0  3689 


RESULTANT  FORCES  ON  STRUCTURE 

LATERAL  FORCE  UPLIFT  FORCE  LATERAL  FORCE 

HEADUATER  SIDE  TAILUATER  SIDE 

(LBS)  (LBS)  (LBS) 


5836  4  9363  0  1779  6 


DO  YOU  UANT  TO  PLOT  UATER  PRESSURES?  YES  OR  NO 
-V 


Figure  C4a.  Program  output  for  Example  2 
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J  i  i  Hi  I  1 


Example  '3,  Dam  Wit-h  One  fr-ft  Sheetpile  at  Heel 

6.  As  shown  in  Figure  (:2c,  Example  3  is  a  dam  with  one  6-ft  sheetpile 

at  heel-  The  method  of  fragments  solution  is  shown  in  Figures  C5a  and  C5b. 

A  plot  of  the  uplift  pressures  for  each  method  is  shown  in  Figure  C5c.  The 

uplift  comparison  is  shown  in  Table  C5a.  The  total  flow  and  exit  gradient 

are  compared  in  Table  C5b. 

COHFPCI  rtoo  -  METHOD  OF  FRAGMENTS 

TINE-  II '84 ■  7  DATE*  1/11/13 


TITLE  -  Daa  «1U  »n  •'  SHaatpila 

0  •  SI  SOS?  (FT»2/DAV> 

K  •  H  (FT/DAV) 

0/K  •  S  66  (FT) 

TOTAL  HEAD  LOSS  •  It  M  (FT) 

FRAG  FRAG 
HO.  TYPE 


L 

(FT) 


A 

(FT) 


■ 

(FT) 


T 

(FT) 


SI 

(FT) 


2 

3 


26. M 


30  00 

3t.M 


EM 
6  M 


UARHIHC  -  THERE  IS  HO  EH1EDHENT  ON  THE  TAILUATER  SIDE 
EXIT  GRADIENT  IS  INFINITE  AND  PIPING  IMV  OCCUR 


S2 

(FT) 


FORN 

factor 


•  62 
1  15 


HEAD 

LOSS 

(FT) 


SI 

48 


RESULTANT  FORCES  ON  STRUCTURE 


LATERAL  FORCE 
HEADUATER  SIDE 
(LtS) 


UPLIFT  FORCE 

(LBS) 


LATERAL  FORCE 
TAILUATER  SIDE 
(LIS) 


7336.8 


3117  6 


3274.6 


DO  YOU  WANT  TO  PLOT  WATER  PRESSURES?  VES  OR  NO 


Figure  C5a.  Program  output  for  F-xample  3 


LENGTH  Of  PRESSURE  PRISM 

4.  a.  ia. 

14  . 


10. 


Fxample  4,  ban:  Kith  Two  h  -  f  t 


Shee t  p i  i es 

6.  As  shown  in  Figure  Cdd ,  Fxample  4  is  a  dam  with  two  f>— I  t  ■  hce  t  p  i  1  •.  ■- . 
The  method  of  fragments  solution  is  shown  in  Figures  fhn  and  C6b.  A  plot  oi 
the  uplift  pressures  for  each  method  is  shown  in  Figure  CCc.  The  u p 1  i : t 
comparison  is  shown  in  Table  (ha.  The  total  f  low  and  exit  gradient  art- 
compared  in  Table  C6b. 


CONFINED  FLOU  -  METHOD  OF  FRAGMENTS 

TIME  1?  10  35  DATE  8/  8/83  £ 

TITLE  -  Daa  with  two  6'  Shoetpllos 
0  -  73  7761  (FTIF3/DAV ) 

K  •  14  4000  (FT/DAV) 

Q/K  •  S  12  (FT)  _ 

TOTAL  HEAD  LOSS  -  10  08  (FT) 


FRAG 

NO 

FRAG 

TYPE 

L 

(FT) 

A 

(FT) 

B 

(FT) 

T 

(FT) 

SI 

(FT) 

S2 

(FT) 

FORM 

FACTOR 

HEAD 

LOSS 
(FT  ) 

l 

2 

30  80 

6  00 

0  62 

3  17 

2 

5 

20  00 

30  00 

6  00 

0  71 

3  65 

3 

2 

30  00 

6  00 

0  62 

3  1? 

-  _ 

D(IT  GRADIENT  •  0  3339 


RESULTANT  FORCES  ON  STRUCTURE 


LATERAL  FORCE  UPLIFT  FORCE 

HEADUATER  SIDE 

(LBS)  (LBS) 


LATERAL  FORCE 

TAILUATER  SIDE  • 

(LBS) 


7393  1 


6240  0 


1717  3 


DO  VOU  WANT  TO  PLOT  WATER  PRESSURES?  YES  OR  NO 
*Y 


Figure  CAn.  Program  output  for  F.xample  4 
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Example  5,  Dam  With  One  15-ft  Sheetpile  at  Toe 


7.  As  shown  in  Figure  C2e,  Example  5  is  a  dam  with  one  15-ft  sheetpile 
at  the  toe.  The  method  of  fragments  solution  is  shown  in  Figures  C7a 
and  C7b.  A  plot  of  the  uplift  pressures  for  each  method  is  shown  in 
Figure  C7c.  The  uplift  comparison  is  shown  in  Table  C7a.  The  total  flow  and 
exit  gradient  are  compared  in  Table  C7b. 


CONFINED  Flow  -  METHOD  OF  FRAGMENTS 

TINE  1?  17' IS  DATE •  By  »y 83 


TITLE  -  Daa  Milk  «••  IS*  Skaatpila 
0  •  68  7944  (FTS*2/DAV) 

C  •  14  4MI  (FT'DAY) 

Q'K  -  4  22  (FT) 


TOTAL  HEAD 

LOSS  • 

18  88  (FT) 

FRAG  FRAG 
HO  TYPE 

L 

(FT) 

A  1 

(FT)  (FT) 

T 

(FT) 

SI 

(FT) 

S2 

(FT) 

FORN 

FACTOR 

HEAD 

LOSS 

(FT) 

1  3 

2  2 

28  88 

38  #8 

38  88 

IS  88 

IS  88 

1  37 

1  88 

S  78 
4  22 

EXIT  GRADIENT  •  *1686 


LATERAL  FORCE 
HEADUATER  SIDE 
(L8S) 


RESULTANT  FORCES  ON  STRUCTURE 
UPLIFT  FORCE 
(LIS) 


LATERAL  FORCE 
TAILUATER  SIDE 
(LIS) 


1S2S8  « 


18419  7 


8995  8 


DO  YOU  WANT  TO  PLOT  WATER  PRESSURES?  YES  OR  NO 
•  V 


Figure  C7a.  Program  output  for  Fxample  5 


me 


,  :jdu'  l  /  a 


Uplift  Comparison  for  Example  5,  1) am_W uj i 
One  15-ft  Sheetpile  at  Toe 


Total  Uplift  Location  of  Resultant 


%  Difference  from 
Finite  Element 


Met  hod 

Force,  lb 

from  Left,  ft 

Uplift  Force 

Resultant  Local 

MOF 

10420.8 

9.341 

8.41 

0.41  1 

CREEP 

9984.0 

9.167 

4.13 

2.30 

FEM 

9580.2 

9 . 380 

Table  C7b 

Tot  a_l  FI  ow  and  Exit  Gradient  Comparison  For  Example  5 , 
Dam  With  One  1  5- ft  Sheetpile  at.  Toe 


Met  hod 
Mi  IF 
FEM 


Total  Flow 
_  ft  V day 

60.8 

65 . 5 


Exit  Grad ient 
0. 1686 
0. 16  16 


%  Difference  from 

Finite  Elemen t _ 

Flow  Exit  Gradient 


Example  6,  Dam  With  One  15-ft  Sheetpile  at  Heel 

8.  As  shown  in  Figure  C2f,  Example  6  is  a  dam  with  one  15-ft  sheetpile 
at  the  heel.  The  method  of  fragments  solution  is  shown  in  Figures  C8a  and 
C8b.  A  plot  of  the  uplift  pressures  for  each  method  is  shown  in  Figure  C8c . 
The  uplift  comparison  is  shown  in  Table  C8a.  The  total  flow  and  exit  gradient 
are  compared  in  Table  C8b. 


CWIMD  now  - 

TIME*  11 <31 '87 


KTNOO  Of  FRACHENTS 


DATE-  1/11/13 


TltlC  -  Baa  Milk  m*  IS'  SkMlplU 
•  •  HUM  IFTttl/MV) 

K  •  14  4MB  (FT'DAY) 

It  ■  4  23  (FT ) 

TOTAL  HIM  LOSS  •  10  M  (FT) 

FANS  FMC 


L  A  1 

I  T 

SI 

nr>  (Ft)  (i 

FT)  (FT)  (1 

FT 

FORI) 

FACTOR 


NARH1NS  -  TICK  IS  NO  CNKMCMT  ON  THE  TAIUIA1SR  SIDE 
EMIT  MAD1EMT  IS  IKINITC  AND  FIFING  NAT  OCCUR. 


LATERAL  FORCE 
KADUATER  SIK 
(LBS) 


RESULTANT  FORCES  OH  STRUCTURE 

ICE  UPLIFT  FORCE  Li 

IK  Ti 

(LBS) 


LATERAL  FORCE 
TAILUATER  SIK 
(LBS) 


175242 


11268.4 


BO  YOU  WANT  TO  PLOT  WATER  PRESSURES?  YES  OR  NO. 


Figure  C8a.  Program  output  for  Example  6 


Table  C8a 


Uplift  Comparison  for  Example  6,  Dam  With 
One  15-ft  Sheetpile  at  Heel 


/  Difference  from 


Total  Uplift  Location  of  Resultant  _  Finite  Element 


Method 

Force,  ibs 

from  Left,  ft 

Uplift  Force 

Resultant  Location 

MOF 

2059.2 

6.667 

33.90 

17.57 

CREEP 

2496.0 

6.667 

14.97 

17.57 

FEM 

2899.8 

7.951 

Table  C8b 

Total  Flow  and  Exit  Gradient  Comparison  for  Example  6, 
Dam  With  One  15-ft  Sheetpile  at  Heel 


Method 

Total  Flow 

f  t-Vday 

Exit  Gradient 

%  Difference  from 
Finite  Element 

Flow  Exit  Gradient 

MOF 

60.8 

* 

7.44  * 

FEM 

65.5 

* 

*  Infinite  exit  gradient  for  case  of  no  embedment  on  downstream  side. 
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Example  7,  Dam  With  Two  15-ft  Sheetpiles 


9.  As  shown  in  Figure  C2g,  Example  7  is  a  dam  with  two  15-ft  sheetpiles 
The  method  of  fragments  solution  is  shown  in  Figures  C9a  and  C9b.  A  plot  of 
the  uplift  pressures  for  each  method  is  shown  in  Figure  C9c .  The  uplift 
comparison  is  shown  in  Table  C9a.  The  total  flow  and  exit  gradient  are 
compared  in  Table  C9b. 


CONFINED  FLOU  -  METHOD  OF  FRAGMENTS 

TIME  17-13  36  DATE-  8/  8/83 


TITLE  -  Daa  with  two  IS'  Shaatpllas 
Q  ■  <7  6S61  (FTMa/DAY) 

t  •  14  4006  (FT/DAV) 

0/K  •  3  31  (FT) 

TOTAL  HEAD  LOSS  •  10  00  (FT) 

FRAG  FRAG  FORM  HEAD 


NO 

TYPE 

L 

(FT) 

A  B 

(FT)  (FT) 

T 

(FT) 

SI 

(FT) 

S2 

(FT) 

FACTOR 

LOSS 

(FT) 

1 

2 

30  00 

IS  00 

1  00 

3  31 

2 

S 

20  00 

30  00 

15  00 

1  02 

3  38 

3 

2 

30  00 

IS  00 

1  00 

3  31 

EXIT  GRADIENT  •  01323 


LATERAL  FORCE 
HEADUATER  SIDE 
(LBS) 


RESULTANT  FORCES  ON  STRUCTURE 
UPLIFT  FORCE 
(LBS) 


LATERAL  FORCE 
TAILUATER  side 
(LBS) 


179S1  2 


6240  0 


8S68  8 


DO  YOU  WANT  TO  PLOT  WATER  PRESSURES?  YES  OR  NO 
•  V 


Figure  C9a.  Program  output  for  Example  7 


Figure  C9b.  Plot  of  pressures  for  Example 


Up  lit 

t  Compar 

Ison 

for 

Examp  1  e 

7  ,  Dam 

With 

Two 

15- 

ft  Sheetpiles 

% 

Di  f  f 

erenee  Iron) 

Tot  a  1 

Uplift  ! 

.oral  ion 

of  K 

e  s  u  1 

tant 

F  i  u  i 

!  e  F.  1  ement 

Met  hod 

Force 

,  lbs 

t"  rom 

Left 

,  IL 

Up! 

i  i  t  Fo i 

Tt1 

Kesu  1 1  an  t  hi 

MOF 

6240 

.0 

9. 

54  7 

0.058 

2.71 

CREEP 

6240 

.0 

'). 

1  h  J 

0.058 

6.7/ 

FEM 

624  J 

.  6 

9. 

809 

Table  C9b 

Total  Flow  and  Exit  Gradient  Comparison  for  Example  7, 
Dam  With  Two  IS- ft  Sheetpiles 


Method 


1'ota1  Flow 
ft  Vday 


%  Difference  from 
F i i lit  e  _E  1  em e n t 

Exit  Gradient  Flow  Exit  Gradient 


Ml  >  F 


0. 1322 
0.1270 


4.01 


FLM 


47.7 

51.2 


7.08 
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